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SECTION  I 
INTRODUCTION 

Electrode  boundary  regions  play  an  important  role  in 
high  energy  density  discharges  presently  used  for  electrical 
lasers.  Discharge  media  of  particular  interest  are  molecular 
gases  in  glow-discharge  generated  plasmas  as  exemplified  by 
the  C02  laser.  A  quantitative  description  of  the  electrode 
regions  is  required  in  order  to  establish  energy  loadings  and 
geometric  scaling  in  such  discharges  for  various  practical 
devices  presently  under  development. 

The  electrode  regions  correspond  to  those  portions  of  the 
plasma  that  are  affected  by  the  presence  of  the  electrodes.1 
Under  certain  conditions,  the  relevant  characteristic  lengths 
of  the  problem  are  independent  of  and  smaller  than  typical 
interelectrode  dimensions.  The  variables,  therefore,  display 
a  "boundary  layer"  behavior.  This  behavior  is  already  well 
understood  for  any  fluid  dynamic  lengths  that  may  be  of 
importance. 2  In  both  stationary  and  moving  plasmas  of  moder¬ 
ately  high  density,  the  electrode  regions  are  of  the  boundary 
layer  type  and  the  electrode  contributions  may  be  described 
with  the  continuum  or  classical  fluid  dynamical  equations. 

In  addition  to  diffusion3,  these  equations  must  describe  the 
electrical  conduction  term,  Poisson's  equation,  and  suitable 
energy  relations.  To  the  extent  that  the  presence  of  non¬ 
equilibrium  effects  may  be  well  represented,  the  continuum 
approach  may  be  both  simpler  than  present  kinetic  theory 
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approaches4 » 1  and  accurate.  Any  classical  approach  brings 
together  a  vast  body  of  literature  unmatched  by  kinetic 
theory,  even  if  the  latter  is  more  appropriate  for  non¬ 
equilibrium  situations.  The  method  in  our  work  corresponds 
to  a  hybrid  approach  to  the  problem. 

Electrode  voltage  drops  are  significant  in  glow  dis¬ 
charges  because  in  addition  to  the  necessary  allowance 
required  in  the  external  power  source,  "one  can  never  impose 
a  uniform  electric  field  on  a  plasma  by  means  of  collecting 
electrodes".6  Thus,  in  the  ionizer/sustainer  discharge  mode 
of  pumping  a  laser,  the  proper  E/N  value  is  obtained  after 
the  electrode  regions  have  been  accounted  for.  Moreover,  the 
values  of  E/N  are  higher  in  the  electrode  regions  than  in  the 
plasma  core  (i.e.,  in  the  positive  column  or  undisturbed 
plasma) ,  and  the  discharge  stability  to  glow  collapse  may  be 
governed  by  what  E/N  values  can  be  tolerated  in  the  electrode 
regions  without  exceeding  some  ionization  rate  critical  value 
which  will  lead  to  instabilities  and  arcing.8  Another  phe¬ 
nomenon  that  affects  both  plasma  homogeneity  and  stability  is 
the  presence  of  constriction  or  spots  at  the  electrodes9 
surfaces.  In  the  so-called  vacuum  arc,10*11  anode  spots 
appear  under  certain  operating  conditions  and  are  of  some 
practical  importance;  anode  spots  have  also  been  observed12 
in  high  density  plasmas.  Since  intense  vacuum  arcs  produce  a 
dense  metal  vapor  at  the  surface,  they  too  are  described  by 
the  continuum  equations  and  ultimately  it  should  be  possible 
to  reconcile  results  from  the  arc  literature  to  our  work. 
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The  intense  cathode  spot  is  the  conventional  mode  of  operation 
of  many  arcs;  however,  we  shall  be  only  concerned  with  describ¬ 
ing  the  "cold  cathode"  because  of  its  relative  simplicity  and 
its  closer  kinship  with  anode  behavior. 

Thus,  the  imposition  of  a  homogeneous,  precisely  defined 
value  of  F/N  in  an  E-beam  generated  plasma  is  not  without 
challenge.  The  E-beam  generated  ionization14  is  usually  less 
homogeneous  than  the  sustainer  field  distribution  so  that  E/N 
and  ne  may  vary  in  space  and  time  throughout  the  discharge. 

For  CW  devices,  we  must  add  fluctuations  in  density  due  to 
turbulence  and  to  heating,  and  in  pulsed  devices  fluctuations 
in  density  due  to  shock  waves.  Fortunately,  in  the  Hg-Ng-CC^ 
laser,  the  optimum  E/N  for  laser  efficiency  is  rather  broad. 

1 .  Approach 

Our  description  of  the  anode  and  the  (cold)  cathode  is  based 
on  the  existence  of  both  a  continuum  and  local  thermodynamic 
equilibrium  (LTE)13.  As  such,  except  for  photons,  particles  are 
locally  at  equilibrium  and  the  Boltzmann  relation  between  energy 
states  may  be  used.  Moreover,  we  sometimes  assume  that  both 
electrons  and  heavy  particles  have  Maxwellian  velocity  distri¬ 
butions  but  at  different  temperatures.  Whenever  electron  colli¬ 
sions  are  the  primary  source  of  transitions,  then  the  electron 
temperature  appears  in  the  corresponding  Boltzmann  relation. 

Thus,  the  non-equilibrium  situation  represented  is  one  where 
the  electrons  have  a  Maxwellian  distribution  at  a  corresponding 
temperature  greater  than  the  translational  or  vibrational  temper¬ 
atures  of  the  lasing  molecules.  The  term  "two-temperature  plasma" 
is  often  used  to  describe  the  above  situation13. 
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At  this  stage,  there  are  a  number  of  simplifications  incor¬ 
porated  in  our  description.  First,  the  effects  of  a  magnetic 
field  are  neglected  since  no  external  magnetic  fields  are 
present  and  since  in  the  unconstricted  mode  the  internal  or 
self -generated  magnetic  fields  is  expected  to  have  no  effect  on 
the  discharge.  Second,  we  consider  the  overall  gas  temperature 
to  remain  constant;  this  is  an  adequate  assumption  for  convec- 
tively  cooled,  high-volume  discharges.  In  such  a  thermally 
stable  system,  however,  transitions  from  a  glow  to  arc  filaments 
cannot  be  represented.  Third  and  last,  the  molecular  gas 
described  is  pure  nitrogen  rather  than  the  required  laser  mixture, 
and  the  source  of  ionization  through  an  E-beam  or  other  source11* 
is  external  to  the  sustainer  power;  the  only  function  of  the  elec¬ 
trodes  is  that  the  imposing  an  E/N  to  pump  the  laser  medium. 

The  equations  that  describe  the  electrode  voltage  drops 
comprise  species  continuity,  the  electron  energy  equation,  and 
Poisson's  equation15’16.  The  small  electron  mass  yields  important 
consequences,  two  of  which  are  the  ability  to  neglect  convection 
in  the  electron  continuity  equation  and  the  strong  dependence  of 
the  electron  temperature  on  E/N.  In  this  report,  we  de¬ 
couple  the  presence  of  the  E-beam  by  assuming  that  the  sustainer 
operates  only  during  the  afterglow  region. 

In  glow  discharges,  the  electron  energy  distribution  is  known 
to  be  non-l'axwell  iar. 1  7  - 1  9  and  suitable  account  must  be  made  of  this 
important  fact.  In  order  to  allow  for  non-Maxwel 1 ian  distribu¬ 
tions  and  still  retain  the  use  of  the  concept  of  temperature. 
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empirical  information  is  used  in  the  representation  of  the  coeffi¬ 
cients20;  thus,  the  electron  temperature,  the  ionization  and 
recombination  coefficients,  the  diffusion  coefficients,  etc., 
are  given  as  parametrized  functions  of  E/M  from  experimental 
measurements . 

In  the  usual  boundary  layer  problem  in  fluid  mechanics  a 
two-dimensional  flat  plate  is  assumed2.  It  is  shown  herein 
that  similarly  a  two-dimensional  Cartesian  description  is  the 
minimum  suitable  description  for  our  equations  in  the  limit  of 
low  currents.  We  further  argue  that  convection  can  only  be 
significant  in  the  plasma  portion  of  the  ambipolar  regions. 

Note  that  compared  to  the  sheath,  the  ambipolar  regions  contribute 
only  a  minor  fraction  of  the  voltage  drop.  The  sheath  and  ambi¬ 
polar  regions,  however,  are  important  in  questions  of  stability 
of  the  discharge  and  must  be  included  in  the  problem  formation. 

We  model  a  quiescent  plasma  in  the  immediate  neighbourhood 
of  the  electrode  surfaces.  In  an  attempt  to  describe  CW  devices, 
we  have  included  in  this  report  a  discussion  of  possible  effects 
of  convection.  In  actuality, no  calculations  of  such  effects  are 
performed.  Now  in  pulsed  lasers,  where  pulse  lengths  of  the 
order  of  50  usee  are  of  interest,  convection  will  play  no  direct 
role.  Of  considerable  importance,  however,  is  the  fact  that 
both  the  sheath  and  ambipolar  regions  can  establish  themselves 
within  such  short  times;  this  possibility  is  a  direct  result  of 
the  boundary  layer  nature  of  the  problem. 
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2. 


Parameters  of  Discharge  Gas 

Throughout  this  work  we  deal  exclusively  with  nitrogen  dis¬ 


charges  under  the  following  coefficients: 

Gas  Density:  1  Amagat  or  3xl025  m~3 

Gas  Type:  molecular  (nitrogen) 

Gas  Temperature:  273  °K 

Electric  Field:  1.5  to  15x10s  V/m  (E/N  =  0.5  to  5xlO~20 

Vm2 ) 

Electron  Density:  1017  to  101  9  m~ 3 

Electrode  Voltag  Drop:  35V 

Since  cross-section  data  for  electron  in  nitrogen  discharges 
are  readily  available18,  the  electron  temperature  and  diffusion 
coefficient  are  calculable.  Also,  n*/N  where  n*  is  a  stationary 
state  (see  Section  II-2)  and  N  the  gas  number  density  can  be 
computed.  These  parameters  are  shown  in  Figure  1  as  a  function 
of  E/N.  The  ions  are  assumed  to  follow  the  gas  temperature .  An 
electrode  voltage  drop  of  about  35  volts  may  be  assumed  since  it 
represents  the  effective  ionization  energy8  for  N9  and,  therefore, 
the  anode  fall.  However,  results  shown  in  Sec.  IV  are  presently 
restricted  to  10  volts  because  of  computat ional  difficulties  at 
the  higher  voltages. 

"The  cathode  is  more  complicated  to  describe  because  of  the 
requirement  of  electron  emission  from  the  surface.  The  cathode 
fall  in  cold  cathodes  is  of  the  order  of  250  volts  and  in  thermionic 
cathodes  considerably  less.  A  proper  representation  of  the  cathode 
will  require  a  more  efficient  computational  scheme  with  suitable 
boundary  conditions. 
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SECTION  II 


PROBLEM  FORMULATION 


1 .  Equations 

The  species  equations  for  singly  ionized  atoms  without 
any  magnetic  fields  may  be  written  as  are  commonly  in  mass 
transfer  analysis1 5 , 1 p  >  2 1 


ion.,:  pu-VC.-y-  [PD.(-  ^  +  ^  Vp.) 


=  W  . 


(1) 


me 

Electrons  :  - - o  -  V- 

e  1 


eC  F.  C 


=  w  ( 2  ) 
e  v 


! 

I 


r 


where 


u  = 


C. 


i  ,e 


w . 


1  ,  o 


D. 


l  .  e 
T 


T  = 

o 


Pi,e  = 
e  = 


me 
k  = 


q  = 


overall  gas  density 

overall  gas  velocity  field 

species  mass  fraction  L.C.  =  1  (pC.=m.n.) 

3  3  J  J  J 

net  source  term  £  .  \V .  =  0 

J  J 

diffusion  coefficient 

electron  temperature 

overall  gas  temperature 

species  partial  pressure  (p^.  =  n^kT^) 

charge  of  the  electron 

mass  of  the  electron 

Boltzmann’s  constant 

electric  field 

rate  of  t'iermal  iznt  ion  of  fast  electrons 
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Futhermore , 


\V’i  =  we  for  single  ionization  in  the  gas  (3) 

and  E  =  -V<t>  where  <t>  is  the  electric  potential  (4) 


In  our  notation,  n  represents  species  number  density 
and  N  the  overall  number  density.  Also,  since  T  is  a 
given  function  of  E/N,  there  is  no  further  need  for  an 
electron  energy  equation.  The  species  mass  fraction,  C, 
is  in  itself  dimensionless  but  it  will  prove  useful  to  define 
a  new  ratio  for  this  variable  (see  Section  III).  Equations  1 
and  2  are  written  for  steady  flow. 

The  only  equation  needed  to  complete  the  set  now  is 
Poisson's  equation, 


(C. - — C  ) 

eomi  1  me  6 


(5) 


where  =  permittivity  of  free  space. 

From  this  set  of  equations,  the  sheath  and  ambipolar 
regions  evolve  in  a  self-consistent  way,  obviating  the 
requirement  to  match  boundary  conditions  between  the  regions. 
2.  Characteristic  Lengths 

It  turns  out  that  all  the  dependent  variables  in  this 
problem,  namely,  Ci ,  C  ,  and  V<J> ,  can  be  of  the  boundary- 
layer  type.  As  such,  there  is  an  individual  characteristic 
length  over  which  the  magnitude  of  these  variables  changes 
from  their  value  at  the  electrode  to  their  value  in  the 
undisturbed  plasma.  Of  course,  the  fluid  dynamic  boundary 
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layer  is  established  by  non-electrical  considerations  and  we 


shall  assume  that  it  ranges  from  one  or  a  few  millimeters  to 
a  few  centimeters  depending  on  the  Reynolds  number  of  the  flow. 

The  concentration  boundary  layers  are  not  easily  surmised 
because  of  inflections  at  the  sheath  and  because,  as  will  be 
pointed  out ,  their  boundary  layer  nature  depends  strongly  on 
the  degree  of  reaction  in  the  plasma,  i.e.,  on  wg  . 

The  extent  of  the  sheath  is  one  of  the  most  important 
characteristic  lengths  in  this  work  because  it  is  within  the 
sheath  that  most  of  the  voltage  drop  is  expected  to  occur  (for 
a  short  discharge).  Fortunately,  the  sheath  length  can  be 
estimated  rather  easily  from  Poisson's  or  in  this  case  Gauss' 
equat ion , 


V.E  =  f-(n.-n  ) 

Eo 


(6) 


where  eQ  is  the  permittivity  of  free  space. 


Let 


(E-E  ) a 

E  =  ttk - r — i—  and  n  . 

!  4>  -<!>  I  i ,  e 


n  . 

1 ,  e 


(7) 


where  Ag  =  sheath  characteristic  length;  the  subscript  "s” 
indicates  the  value  at  the  electrode  (anode  or  cathode)  and 
the  value  at  the  undisturbed  plasma.  Equation  7  then 
becomes 


V  •  E  = 


en»V 


(n . -n  ) 
i  e 


(8) 
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In  the  equation  above,  A  is  used  to  non-dimensionalize 

s 

the  del  operator  in  all  directions.  That  is,  independent  of 
the  dimensionality  of  the  problem,  the  sheath  is  effective 
over  the  characteristic  dimension  of  Ag  .  Now  we  can  esti¬ 
mate  the  size  of  the  sheath  by  setting  the  square  brackets  in 
Equation  8  to  be  of  order  one:  z,  since  all  other  terms  have 
been  made  of  order  one, 


where  A^  =  the  Debye  length  in  the  undisturbed  plasma. 

We  now  establish  the  specific  values  anticipated  for 
Ag  .  Take  for  example 

n  =  1017  to  1019m~3 

00 

=  35  VOltS 

Then,  Ag  =1.4x  10-4  toi  qx  10-5m 

It  is  interesting  to  note  that  a  change  of  I I  by 
a  factor  of  10  would  change  Ag  by  a  factor  of  3  so  that 
we  may  consider  the  estimate  reasonable  for  the  above-quoted 
electron  number  densities  in  the  undisturbed  plasma.  Now,  if 
the  fluid  dynamic  boundary  layer  for  a  density  of  one  Amagat 
is  at  least  an  order  of  magnitude  greater  than  the  sheath, 
then  we  may  conclude  that  convection  will  have  a  small  if  not 
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negligible  effect  in  the  sheath.  Note  that  the  sheath,  being 
always  adjacent  to  the  surface,  resides  in  the  region  where 
the  velocity  drops  to  zero;  therefore,  it  is  entirely 
appropriate  to  assume  that  convection  plays  no  significant 
role  in  the  sheath. 

Next  in  importance  is  the  length  of  the  ambipolar  region, 
i.e.,  the  transition  region  between  the  sheath  and  the 
undisturbed  plasma.  This  region  is  neutral  and  can  perhaps 
span  the  boundary  layer  so  that,  in  CIV  devices,  convection 
would  likely  be  present.  We  do,  however,  expect  that  convec¬ 
tion  may  be  neglected  without  affecting  the  resulting  voltages 
appreciably.  In  Equation  1,  for  example,  we  know23  that  in 
the  ambipolar  region  the  last  term  of  the  left  hand  side  is 
small  compared  to  electric  conduction;  therefore,  we  need  only 
compare  convection  with  conduction,  or  the  fluid  velocity  with 
the  drift  velocity.  The  transverse  velocity  component  in  the 
boundary  layer  problem  is2 

v  ~  U*,  Re““  ~  0.1-lm/s 

assuming  the  Reynolds  number  to  be  Rq  ~  ] 0 4  -  10c  and  the 
free  stream  velocity  to  be  Uto  ~  10  -  100  m/s  .  Now  the 
drift  velocity  may  be  approximated  by 

v„  ~  y.E  ~  10  -  100  m/s 

for  E^  1'  10"  -  1 0 6  V/m  and  10-4  m2/s-volt.  Clearly, 

in  the  presence  of  a  sufficiently  strong  interelectrode  field, 
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the  contribution  of  ion  convection  due  to  a  cross  flow  should 
be  negligible.  The  geometrical  orientation  of  the  discharge 
with  respect  to  the  flow  is  an  important  factor  here.9 

The  concentration  profiles  behave  in  a  more  complicated 
manner  than  either  the  voltage  or  the  electric  field.  While 
the  magnitude  of  and  np  will  change  appreciably  within 

the  sheath,  fractional  analysis  is  risky  because  inflections 
are  present  in  the  profiles.  In  this  section  and  in  Appendix 
A,  we  give  a  discussion  pertaining  to  the  stability  of  the 
ambipolar  region  and  its  boundary  layer  nature.  Since 
(see  Equation  A6  and  Figure  1),  the  characteristic  length 
for  ambipolar  diffusion,  turns  out  to  be  comparable  to  the 
sheath  length,  ambipolar  diffusion  to  the  walls  can  establish 
itself  during  periods  of  the  order  of  50  u  sec. 

In  general,  the  ambipolar  diffusion  equation  can  be 
written  as 

-  D  V2n  =  ft  (10) 

where  Da  =  ambipolar  diffusion  coefficient 

If  we  may  assume  a  form  for  the  net  production  rate, 
ne  ,  some  important  conclusions  may  be  drawn  about  the  sta¬ 
bility  of  the  ambipolar  region.  Let 

n  =  v  .  n  -  an 2  +  '!'  (11) 

e  i 

where  =  ionization  rate  coefficient,  sustainer 

a  =  three-body  neutral  recombina¬ 
tion  rate  coefficient 

t  =  ionization  rate  coefficient  E  beam 
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In  a  nitrogen  discharge,  no  electron-neutral  attachment  is 
present.  (Presumably,24  attachment  instabilities  can  be 
significant  in  lasers,  but  are  not  described  here.) 

If  we  look  now  at  the  steady  form  of  the  equation  we  have 


V2n  =-g£(n*_n)  -  ¥/Da 


where 


i  *  = 


v .  /  a  or  n* /v  -■ 

1'  'a 


v./N 


(12) 


The  parameter  n*  above  is  governed  solely  by  the  sustainer 
discharge,  i.e. ,  by  E/N  as  shown  in  Figure  1.  It  is  reasonable 
to  assume  that  E/N  is  designed  to  be  fairly  homogeneous  in  the 
interelectrode  space  but  that  the  E-beam  produced  charge  con¬ 
centration  is  not14.  Now  at  the  outer  edge  of  the  ambipolar 
region  (at  the  undisturbed  plasma),  n(x,y)  =  n^  ,  where  nTO  is 
assumed  to  be  governed  exclusively  by  the  E-beam  conditions. 

If  we  now  assume  that  we  operate  in  the  afterglow  region. 

anOT 

=~j5  (n*-nTO)  (13) 

a 


From  the  properties  of 
the  following: 

i)  for  n*  >  nOT,  V2nOT 
possible;  n  can 
ii)  for  n*  <  n^  ,  V2n 
not  possible;  n 


the  Laplacian  operator,  we  can  infer 

<  0  and  a  steady  solution  is 
reach  a  maximum. 

>  0  and  the  steady  solution  is 
cannot  be  a  maximum. 


V 


.*• 
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iii )  for  n*  =  nOT  ,  n  can  neither  be  a  maximum  nor  a 
minimum  inside  the  domain  and,  indeed,  it  must  be 
a  constant . 

The  reason  that  the  second  case  is  impossible  is  that  in 
the  ambipolar  region  V2n  <  0  and  n  seeks  a  maximum  toward 
the  centerline.  If,  however,  n*  <  n^  then  this  trend  must 

reverse  itself  somewhere  and  a  minimum  must  exist  within  the 
domain.  Such  a  situation  corresponds  to  a  physically  unsteady 
condition  and  the  full  form  of  Equation  10  must  be  investigated. 

The  criterion  for  stability  appears  to  be  n*  >  nm  ,  with 
the  equal  sign  as  the  marginally  stable  case.  n^  is 
considered  to  be  a  steady  distribution  attributed  solely  to 
the  E-beam;  this  simplification  of  the  physics  models  the 
ionizer/sustainer  as  if  the  plasma  is  externally  generated 
and  a  steady  distribution  (  i^)  is  independent  of  E/N  and 
reasonably  constant  for  the  pulse  duration. 

In  pulsed  lasers,  diffusion  may  only  be  established 
within  the  boundary  layers.  This,  however,  is  in  itself 
significant  because  the  electric  field  reaches  a  maximum 
within  the  electrode  regions.  But  the  problem  cannot  be 
simplified  here  and  one  must  work  with  the  entire  set  of 
equat ions . 

The  use  of  the  ambipolar  diffusion  equation,  i.e., 

Equation  10,  is  the  traditional  approach  to  the  problem2'. 

This,  however,  rather  complicates  the  calculation  of  the 
electric  field  as  it  does  in  the  Cchottky  solution. 2 ^ 

This  complication  does  not  seem  to  appear  when  the  full 
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formulation  of  the  problem  including  the  sheath,  is  utilized. 
In  order  to  get  some  intuition  into  the  nature  of  the 
solutions,  it  is  worthwhile  to  attempt  some  sort  of  simpli¬ 
fication.  In  Appendix  A,  the  ion  flux  equation  in  the 
ambipolar  region  is  investigated;  it  is  assumed  that  this 
equation  is  sufficient  to  describe  the  concentration  profile 
given  an  electric  field  and  the  form  of  ng  .  We  further 
assume  that  a  one-dimensional  description  is  possible  and, 
indeed,  desirable  since  it  yields  a  constant  electric  field 
in  this  region.  Results  indicate  that  for  E/N  values  of 
interest  the  ambipolar  charge  concentration  is  indeed  of 
the  boundary  layer  type. 

3.  Dimensionality  of  the  Problem 

In  the  limit  of  low  currents,  the  elevation  of  electron 
temperature  may  be  neglected  together  with  ionization  from 
the  sustainer  discharge,  and  the  quiescent  plasma  governing 
equations  may  be  written  in  a  simpler,  equivalent  form.  We 
shall  further  examine  the  one  dimensional  form  because,  at 
first  glance,  the  smallness  of  the  sheath  compared  to  a 
typical  electrode  dimension  suggests  that  the  one-dimensional 
Cartesian  description  might  indeed  be  adequate.  As  before, 
y  is  the  interelectrode  coordinate. 


e2D 


3  i  kT 


o 


dn  . 

-  n.E  -  eD. 

1  i  dv 


e2D 


Je  kT 


o 


dn 

e  e 

—  n  E  +  eDf- 
e  r- 


dv 
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(16) 


(17) 


Mow  in  one-dimensional  flow  with  no  ionization  or  recombina¬ 


tion,  j  ^  and  j  are  individually  constant  throughout  the 
interelectrode  space. 

Subtracting  Equation  17  from  16  we  obtain 


Now  let  us  look  at  the  sign  of  the  terms  in  the  above 
equation  for  three  regions,  namely,  (a)  in  the  anode  sheath, 
(b)  in  the  undisturbed  plasma,  and  (c)  in  the  cathode  sheath. 


Region 

r-=i—  ( n  .  -n  ) 
kT  v  l  QJ 
o 

_~ay(ni  +  ne} 

K 

(a) 

negative 

negat ive 

negat ive 

(b) 

0 

0 

0 

(c) 

posit ive 

posit ive 

posit ive 

Clearly , 

i  1  J  i  >  J  g  » 

,  and  are  constants 

then  Equation  19  is  invalid.  What  is  needed  is  for  j  ^  and 
j  to  decrease  from  the  anode  to  the  plasma  and  to  increase 
from  the  plasma  to  the  cathode  thereby  producing  current 
concentrations  at  the  electrodes12.  It  is  therefore  obvious 
that  for  the  frozen-property  flow  of  current  a  one-dimensional 
formulation  is  inadequate. 
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In  order  to  see  if  ionization  is  a  requisite  for  a  one¬ 
dimensional  solution  we  use  the  conservation  equations  in 
the  form  of  Equation  19  and  we  get 


d_ 

dy 


ne  (D.  +  D  ^ 
1  e 


(20) 


Bui  since  D  >>  D. 

e  l 


we  may  simplify, 

dK~~  % 

dy  D. 

l 


(21) 


Figure  2  shows  the  behavior  of  K  as  defined  in 
Equation  19.  As  can  be  seen  from  the  sketch,  the  term  has 
either  a  positive  or  zero  slope,  so  we  may  conclude  that 

w  or  n  >  0  (a;  anode  region 

v 

=  0  (b)  undisturbed  plasma 

>  0  (c)  cathode  region 

Of  course,  this  means  that  ionization  must  exceed 
recombination  by  an  appropriate  amount  in  order  to  satisfy 
Equations  16-18  and  10,  11.  As  mentioned  previously,  we 
shall  assume  that  ionization  is  due  to  electron  impact  and 
that  recombination  is  of  the  two  or  three-body  type  as 
appropriate  for  moderately  high  pressure  discharges. 

The  material  presented  in  this  section  points  out  some 
interesting  facts.  These  may  be  summarized  by  stating  that 
the  geometry  of  the  current  flow  is  not  necessarily  imposed 
by  the  electrode  geometry.  Depending  on  the  level  of  the 
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Figure  2.  Behavior  of  K-  at  the  Anode,  Undisturbed 

Plasma  and  Cathode  Regions 
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current,  the  plasma  constituents,  etc.,  we  may  have  a  spot 
mode  or  a  glow  mode  (one-dimensional  disregarding  end  effects) 
at  one  or  both  of  the  electrodes12.  Moreover,  the  sheath 
and  ambipolar  regions  may  grow  with  increasing  current 
making  the  solution  of  the  problem  a  challenging  one.  It  is 
clear  that  one  does  not  a  priori  specify  the  dimensionality 
and  size  of  the  domain  for  calculation  but  that  one  has  to 
make  certain  that  the  description  will  be  sufficiently  un¬ 
restricted  to  permit  a  solution  of  the  problem. 

There  are  other  instances  where  a  given  formulation  is 
tractable  only  for  certain  geometries  or  shapes  such  as  the 
flow  of  a  uniform,  incompressible  viscous  fluid  due  to  a 
moving  body  at  small  Reynolds  numbers  (Stoke ' s  flow)27. 
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SECTION  III 
PROBLEM  SOLUTION 

1.  Working  Equations 

It  has  been  shown  above  that  the  depending  on  the  level 
of  the  current  and  associated  voltage,  the  electrode  regions 
represent  a  multidimensional  problem.  We  are,  therefore, 
modeling  non-emitting  electrodes  as  depicted  in  Figure  3. 

Here  a  periodic,  2-dimensional,  flat  plate  region  is  seen. 

We  assume  the  coefficients  to  be  either  constant  as  in  the 
case  of  the  ions  and  neutrals} or  dependant  on  E/N.  For  the 
range  of  E/N  considered,  it  is  adequate  to  take  the  electron 
diffusion  coefficient  and  the  electron  temperature  as 1 8  > : 8 


D 

T 


e 


e 


5.5  x  10  2  (m2/s)  (22) 

12 . l&n(E/NxlO+2  0 )  +  38, 5t  o 
38.5  e~ 


Note  that  E/N  is  in  Volts-m2  in  the  above. 

Since  we  are  neglecting  convection  and  assuming  that 
the  sustainer  operates  in  the  afterglow  of  the  E-beam  current, 
the  governing  equations  in  our  model  are  written  as  follows 


A  _  A  A  A  A  A  0  A 

-V-  lRn.V<j>  +  Vn.j  =  ~  n 

l  l  D.  e 

A  A  1 

n  V<j>  aa  n 

"  Vne  "  IT  VC  " 


(24) 

(25) 


V2<f>  = 


Yp(ne  “ 


Hf) 


(26) 
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Figure  3.  Anode  Model  with  Periodic  Current  Const 
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Now,  we  specialize  our  formulation  to  the  flate  plate 
depicted  in  Figure  3.  We  drop  the  carets  (*)  for  simplicity, 
but  it  is  to  be  understood  that  all  variables  have  been 
suitably  non-dimensionalized .  In  particular, 


y6ni(ni-ne)  -  8 


(3Ji  a  /M„ 

v* 5x  3y  9y/  \3x2  ay2  r\Di)  ‘ 
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5  e  1  3  x  3  y  9y  I  B  ne  I  3  x  3  x 

_  a-l  U  9  3  ne  ,  3  6  3  ne\  (d  ,  9  ^e) 

3 y  3y/ "  \3x3x  3^  )\T7~  %" r?"/ 


(27) 


(28) 


23 


(29) 


mu 


1  »w^t* 


The  associated  boundary  conditions  are  shown  in  Figure 
4.  Each  node  is  in  a  periodic  field  whose  side  or  x-dimension 
is  given  by  a  few  sheath  lengths.  This  "packing"  is  arbitrary 
but  yields  a  maximum  total  current  at  the  electrode  for  a  given 
electrode  voltage  drop.  The  latter  is  governed  by  the  individual 
node  current.  Thus,  we  model  a  maximum  "crowding"  condition. 

Any  further  crowding  would  precipitate  the  same  situation  that 
negates  the  one-dimensional  solution.  The  top  or  y-dimension 
is  unbounded,  as  required  by  the  boundary-layer  type  of  behavior. 
Note,  furthermore,  from  Figure  4  that  the  electron  density  at 
the  anode  is  not  zero  as  is  often  done  in  conventional  probe 
analysis.  Such  a  change  is  needed  because  a  finite  node  current 
requires  a  finite  node  charge  density.  The  node  density  is 
found  by  trying  various  values  until  a  match  of  the  current  at 
the  node  and  the  current  at  the  undisturbed  plasma  is  obtained. 

The  continuous  domain  around  one  node  is  modeled  with 
equidistant  grid  'points  on  a  square  mesh  as  shown  in  Figure  4. 

The  finite  difference  equations  approximating  the  two  species 
equations  and  Poisson's  equation  are  solved  simultaneously  in  a 
line-iterative  fashion. 
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B.  C. 

<f>  =  0  Note- 

n  =  ]  9  calculated  separately 

■  from  coupled  equations 

=  1  that  determine  <£,  ng 

and  n.  . 

Figure  4.  Boundary  Conditions  for  Computational  Method 
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2.  Program 

The  electron  and  ion  species  equations  (Equation  27  and 
28)  contain  non-linear  terms  which  present  problems  in  their 
computer  analsyis.  The  Jacobi  method  includes  all  non-linear 
terms  on  the  "right  hand  side",  i.e.,  external  to  the  coeffi¬ 
cient  matrix.  Convergence  to  a  solution  is  possible  if  these 
non-linear  terms  change  slowly  enough  over  each  iteration.  It 
was  found  by  Dolson  in  Reference  29  that  the  Jacobi  method  was 
in  fact  unstable  for  the  present  set  of  equations  and  conditions. 
As  a  consequence  of  the  failure  of  the  Jacobi  method,  a  quasi- 
Jacobi  method  was  used,  in  which  an  estimate  was  computed  for 
each  total  variable  (<ji,  n  ,  n.)  at  each  grid  point  during  each 
iteration.  When  the  product  of  two  variables  was  encountered, 
one  variable  was  treated  as  a  constant  coefficient  for  each 
iteration.  This  means  that  the  non-linear  terms  are  retained 
in  the  coefficient  matrix.  The  "constant"  coefficients  are 
updated  after  every  iteration,  thus  changing  the  coefficient 
matrix.  The  conventional  Jacobi  method  was  found  to  converge 
only  for  very  low  values  of  <f>  at  the  electrode,  whereas  the 
quasi-Jacobi  procedure  provided  converged  solutions  for  values 
of  <f>  <  2  volts  .  Thus,  the  Newton-Raphson  method  was  chosen. 

The  Newton-Raphson  method,  Reference  30,  is  presented  by 
way  of  illustration.  Equations  27,  28  and  29  are  the  three 
non-dimensional  equations  for  <J>  ,  n  ,  and  n^  .  The  Poisson 
Equation  can  be  written  at  grid  node  i,j  as: 


ili  +  111  + 

3x2  3y2 


Y(ni-ne) 

P 


(30) 
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The  solution  in  this  case  is  said  to  have  converged  ideally  if 

h  k 

F,  .  .  =  0  or  practically  if  the  k  iteration  value  F„  .  . 

1 .  iJ  lj ij 

<  e  .  Writing  this  equation  in  a  linear  truncated  Taylor 
series  expansion,  and  dropping  i,j  subscripts  for  simplicity: 


F1k+1  =  F^  +  AF1k+1 


where 


AF1  =  An.) 


,k+l  ,k  ,  .  ,k+l  . 

$  =  $  +  A<p  etc. 


3 2  <tik+1  ,  3 2  <f k+1  ,  ,  k+1 

- i  +  +  y  ( n  . 

3x 2  3 v 2 


k+1 , 
"e  1 


+  iiik+  Y(„Y  -  „  k>  *  +  i!«k+1  +  (A  k+i  _  A 

3x 2  3y 2  1  6  3x 2  3y 2  1  6 


This  equation  must  hold  at  each  of  the  ixj  grid  points  of  the 
finite  difference  mesh. 

A  typical  non-linear  term: 

V^-Vn  =  &  -5-^  +  |i  -5-^  (33; 

e  3x  3x  3y  3v 

The  expanded  iteration  form  would  look  like,  (x-terms  on lv): 


k  a/v„  k  +  l 
~  ,  3 An 

3<j>  _ e 

3x  3x 


k+1 

11 

k+1 

ne 

k 

li 

-  k 

.^e 

k+1 

3n  k 
e  + 

dX 

3x 

3x 

3x 

3x 

3x 

k  +  1 

.  .  k+1 

3An 

e 

.  3A4> 

+  3x 

3x 

I 


AF  is  considered  linear  in  A<t> ,  Ang  and  An^  and  so  the  last 

term  in  equation  34,  is  neglected,  and  the  products  of  the 

t  h 

solutions  are  known,  as  are  the  k1  coefficients  of  the 

.  . , k+1  .  .  k+1 

unknowns,  A<f>  and  Ang 

Thus  the  final  matrix  has  the  form: 


ktk  Coefficients 
Matrix 


—  — 

ek, 

'k+1  ' 

A  <*■ 

F1 

k 

ne  > 

= 

F2 

n.k 

X 

F3 

(35) 


where  the  F^'s  correspond  in  form  to  Equation  30,  for  each  of 
Equations  27,  28  and  29. 

Thus  the  solution  procedure  consists  of  evaluating  the 
right  hand  side  based  upon  the  kth  solution  and  solving  the 
system  of  equations  for  A<})k+\  Angk+1  and  An^k+1  which 
are  used  to  update  their  respective  grid  point  values. 
Appendix  C  shows  the  program  listing  with  additional  detail. 
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SECTION  IV 


RESULTS 

Computer  solutions  were  obtained  for  the  set  of  conditions 
presented  in  Table  I.  The  anode  voltage  is  ] 0  volts  in  all 
cases , 

TABLE  1:  COMPUTER  SOLUTIONS  CONDITIONS 


CASE 

”e 

9 

H 

I 

not  coupled 

not  coupled 

2.0 

II 

coupled 

not  coupled 

2 . 0 

III 

coupled 

not  coupled 

0.5 

IV 

coupled 

coupled 

0.5 

The  results  are  shown  in  Figures  5-8  as  three-dimensional  views 
of  the  potential  field,  temperature  field,  space  charge  density, 
electron  and  ion  densities.  Also  included  are  graphs  of  poten¬ 
tial,  species  and  space  charge  densities  along  a  cut  from  the 
anode  to  the  free  stream  for  each  case. 

Table  I  gives  the  five  case  conditions,  where  repre¬ 
sents  the  "right  hand  side"  of  Equations  (28)  and  (29),  0 

represents  the  temperature  equation  which  is  a  function  of  the 
solution  set  (E/N,  properly)  and  II  is  equal  to  the  size  of 
the  computational  node  spacing,  i.e.,  HAg  is  physical  space. 

Initial  solutions  were  obtained  with  hG  =  0  with  the 
sheath  size  smaller  than  the  grid  spacing.  These  coarse  solu¬ 
tions  were  illustrative  of  the  ambipolar  solution  in  which 
V 2 $  1  0  except  at  the  anode  where  the  space  charge  is  forced 
as  a  "boundary  condition".  The  effect  of  "turning-on"  the 


Figure  5.  Three  dimensional  presentations  of  CASE  I.  (a)  poten¬ 
tial  field,  (b)  temperature  field,  and  (c)  space  charge  density. 
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Figure  5  (continued).  Three  dimensional  presentations  of  CASE  I. 
(d)  back  view  of  electron  density  and  (e)  backview  of  ion  density. 


DISTANCE  (CM) 


0  =  10  volts  ne  not  coupled 

Aj  =  2.35  10" 5  meter  isothermal 

H  =  2.0 


Figure  5  (f).  Two  dimensional  presentation  of  potential 
distance  from  the  anode  to  the  freestream. 


DISTANCE  (CM) 


4>  =  10  volts  ne  not  coupled 

Aj  =  2.35  10"5  meter  isothermal 

H  =  2.0 


Figure  5  (h) ,  Two  dimensional  presentation  of  space  charge 
density  vs.  distance  from  the  anode  to  the  freestream. 
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POTENTIAL  FIELD  TEMPERATURE  FIELD  SPACE  CHARGE  DENSITY 
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POTENT  I RL  (VOLTS) 
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DISTANCE  (CM) 


<6  =  10  volts  ne  coupled 

Hi  =  2.35  10" 5  meter  isothermal 

H  =  2.0 


Figure  6  (f) .  Two  dimensional  presentation  of  potential  vs. 
distance  from  the  anode  to  the  freestream. 
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DISTANCE  [CM) 


<P  =  10  volts  ne  coupled 
A j  =  2.35  10" 5  meter  isothermal 
H  =  2.0 


Figure  6  (g) .  Two  dimensional  presentation  of  electron  and  ion 
densities  vs.  distance  from  the  anode  to  the  freestream. 
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DISTANCE  (CM) 


<t>  =  10  volts  ne  coupled 
A$-  2.35  10"5  meter  isothermal 
H  =  2.0 


Figure  6  (h) .  Two  dimensional  presentation  of  space  charge 
density  vs.  distance  from  the  anode  to  the  freestream. 
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Figure  7  (g) .  Two  dimensional  presentation  of  electron  and  ion 
densities  vs.  distance  from  the  anode  to  the  freestream. 
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Figure  7  (h)  .  Two  dimensional  presentation  of  space  charge 
density  vs.  distance  from  the  anode  to  the  freestream. 
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Figure  8  (f).  Two  dimensional  presentation  of  potential  vs. 
distance  from  the  anode  to  the  free  stream. 
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Figure  8  (g) .  Two  dimensional  presentation  of  electron  and  ion 
densities  vs.  distance  from  the  anode  to  the  freestream. 
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Figure  8  (h) .  Two  dimensional  presentation  of  space  charge 
density  vs.  distance  from  the  anode  to  the  frees tream. 
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production  term  (rt  )  is  to  bring  about  the  boundary  layer 
behavior  of  the  solution  (Case  I I). The  boundary  layer  behavior 
of  the  species  densities  is  easily  illustrated  in  the  compari¬ 
son  of  Figure  5(g)  (Case  I )  to  Figure  6(g)  (Case  II).  Case  I 
and  II  conditions  produced  an  insufficient  electric  field  to 
enable  the  electron  temperature  tt  to  be  greater  than  1. 

Thus  Case  I  and  II  conditions  were  identical  with  the  tempera¬ 
ture  routine  coupled  or  uncoupled. 

A  solution  would  always  be  obtained  as  long  as  the  sheath 
is  not  larger  than  the  grid  size  spacing,  for  ne  =  0  . 

However  reasonable  solutions  could  be  obtained  for  grid  size 
spacing  smaller  than  the  sheath  size  only  for  hp  f  0  . 

Cases  III  and  IV  each  represent  the  sheath  as  2  grid  spacing 
(^  =  0.5).  Case  IV  is  identical  to  Case  III  with  the  addition 
of  the  electron  temperature  effects.  The  temperature  effect 
in  Case  IV  produced  only  a  slight  increase  of  space  charge 
density  and  a  very  small  increase  of  the  electron  field  near 
the  anode. 

In  all  cases  the  current  into  the  inactive  wall  was  3  or 

4  orders  of  magnitude  less  than  the  current  at  the  anode  node. 

The  current  at  the  anode  node  was  approximately  3  times  the 

total  current  at  the  free  stream.  The  anode  current  could  be 

better  matched  to  the  free  stream  current  by  picking  a  suitable 

ne  (anode),  but  this  was  not  pursued  because  previous  results 

of  a  preliminary  nature  indicate  only  minor  changes  over  the 

results  shown.  In  fact  in  cases  I  and  II  with  n  =  0  ,  at  the 

e 

anode,  the  anode  current  still  exceeds  the  current  at  the 


J 
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free  stream  boundary.  The  electric  field  at  the  free  stream 
boundary  is  insufficient  to  offer  a  current  matching  condition. 
A  more  proper  free  stream  boundary  condition  would  be  E  =  E^ 
rather  than  $  =  0  .  The  electron  field  boundary  condition 
created  numerical  instabilities  and  is  recommended  as  a  subject 
for  later  analysis.  Our  approach  thus  far  was  to  increase  <pQ 
large  enough  to  produce  the  proper  Ew  . 


SECTION  V 


SUMMARY  AND  CONCLUSIONS 

The  results  presented  in  this  report  are  based  on  a  two- 
dimensional  description  of  the  current  flow  at  the  anode. 

Diffusion  is  properly  accounted  for  as  is  ionization/recombi¬ 
nation.  No  convection  or  magnetic  field  effects  are  included. 

We  assume  operation  in  the  after-glow  of  the  E-beam  ionization 
so  that  tp  and  q  are  essentially  uncoupled  from  the  solu¬ 
tions.  Current  constricts  at  active  nodes  along  the  surface, 
and  the  sheath  and  ambipolar  region  are  self-generating  in  the 
solution  of  Equations  27-29. 

The  boundary  layer  nature  of  the  sheath  is  clearly  evident 
in  the  results  shown;  the  boundary  layer  nature  of  the  ambi¬ 
polar  region  only  shows  up  when  wg  or  he  are  coupled  into 
the  calculations.  The  magnitude  of  ne  is  highest  within  the 
sheath  and  along  the  walls,  dropping  off  towards  the  undis¬ 
turbed  plasma.  The  space  charge  density  (ne~n^)  peaks  at 
0.5xl018m~3  near  the  active  node. 

About  90%  of  the  voltage  drop  takes  place  within  one  or 
two  sheath  lengths,  this  produces  a  maximum  electric  field  of 
0 . 8x10 5  Y/m  (for  the  10  volts  drop)  at  the  electrode  node.  7110  undis¬ 
turbed  region  electric  field  is  about  0.5xl04  V/m  which  is 
considerably  below  the  value  of  0.27x10s  V/m  required  for 
pumping  the  laser  medium.  But  then  again,  $  is  10  V  and  it 
should  be  35  V.  That  is,  roughly,  an  increase  by  a  factor  is 
54  required  for  E  .  The  temperature  field,  being  a  direct 
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function  of  E/N,  does  not  affect  our  results  simply  because 
E/N  is  yet  too  low.  There  are,  however,  subtle  changes  which 
we  feel  are  going  to  become  important  at  the  higher  voltages. 

It  is  anticipated  that  H  can  be  decreased  further  (up 
to  1/6  of  a  sheath  length),  with  the  coupling  of  ne  ,  and 
that  this  will  allow  smoother  changes  within  the  sheath  and 
more  numerical  stability  at  higher  voltages.  Also,  it  should 
be  possible  to  consider  jointly  the  E-beam  ionization  in  our 
solutions  to  produce  more  realistic  effects. 
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APPENDIX  A 


AMD I POLAR  REGION  —  ONE  DIMENSIONAL  DESCRIPTION 

We  have  seen  that  within  the  sheath,  the  solution  to  the 
conservation  equations  cannot  be  one  dimensional.  While  there 
may  be  some  particular  range  of  ionization/recombination  which 
permits  one-dimensional  solutions,  most  often  we  must  allow 
for  current  concentrations  at  the  electrodes.  However,  away 
from  the  electrodes  -  the  problem  may  become  one  dimensional 
in  the  ambipolar  region  (see  Figures  3,  Section  III.l). 

Since  the  extent  of  the  sheath  is  quite  small  when  compared  to 
typical  interelectrode  dimensions,  it  is  interesting  to  examine 
one-dimensional  solutions  within  and  beyond  the  ambipolar 
region.  Here  we  assume  that  an  initial  concentration,  ng  , 
is  known  at  the  wall-end  of  the  ambipolar  region  and  that  a 
symmetry  plane  exists  for  the  charge  concentration  (usually 
the  mid-plate  region),  see  Figure  A1 . 

Equations  A1 ,  A2 ,  and  A5  are  shown  below  for  the  one¬ 
dimensional,  ambipolar  region.  We  use  number  density  instead 
of  mass  fraction  and  take  S  to  represent  the  E-beam  production 
rate  of  secondary  electrons.  Convection  is  neglected  and  the 
resulting  electric  field  in  this  region  is  constant.  Also, 
all  the  coefficients  turn  out  to  be  constant. 


1 

j 


Anode 


Figure  A1 .  E-Beam  Discharge  Configuration 
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dv  kT 


,  er. 

d  1 


D.  ^ 

i-  dy 


A  +  S 
e 


d  -ene  p 
dy  [  kTen  - 


A  +  S  +  o 
e 


dE/dy  =0  or  E  =  constant 


We  now  study  Equation  A1  in  the  afterglow  region.  Let 
S  and  q  be  zero  so  that 


n  =  v  .  n  -  an 
e  l 


and  combine  with  Equation  A1 


n.  [eF.  ,  .2 

i  dn  d  n 


.  kT  dv 
vi I  o  •  dv 


=  n  -  n  /n’ 


In  order  to  nondiment ional ize  the  above  relation  we  in¬ 
troduce  the  ambipolar  diffusion  characteristic  length  along 
with  suitable  ratios  of  the  variables. 


-  /  D .  K 

D  -  /v  or  L  N  =  X 
i  i  n  yj  v  ^  /  N 


y  =  y  / 1 


e(E/N)Ln?- 


The  resulting  differential  equation  becomes 


n(±l)  =  ns 

It  is  obvious  from  equation  A7  that  if  n'(0)=0  and 
n(0)=l  we  have  truly  boundary-layer  behavior  since  both 

A  A 

n'  and  n"  become  zero  at  the  midplane.  This  represents  a 

mathematical  limit  but  note  that  because  is  very  small 

in  relation  to  the  interelectrode  distance,  n(0)  will  be 

for  all  practical  purposes  equal  to  one,  i.e.,  n~1.0  . 

Figure  1  in  the  text  shows  calculated  values  of  LnN  ,  as 

a  function  of  E/N.  These  compare  reasonably  well  with  values 

25 

quoted  in  the  literature. 

Equation  A7  is  written  below  without  the  carets  and  with 
primes  instead  of  derivatives. 

n"  -  En'  +  (n-n2)  =  0  (AS) 

n  '  (0)  =  0 
n(±£)  =  ng 

We  have  considered  the  above  equation  for  E  5  51  (or 
-1G  2 

E/N  s  10  V-cm  ).  The  solution  of  Equation  A8  was  obtained 
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numerically  and  is  shown  in  Figure  A2  .  We  note  that  this  solu¬ 
tion  can  be,  within  a  very  close  approximation,  easily  obtained 
by  neglecting  n".  We  have 


and 


n'  *  i(n-n2) 
n '  ^  0  when  n  =  1 


(A9) 


The  important  result  is  that  n'  =  0(1/E)  and  n"  ~  0(1/E2)  so 
that  only  for  E  sufficiently  large  may  we  neglect  diffusion  in 
the  ambi polar  region. 

Solving  Equation  A9  we  obtain: 


n 


1  +  (l/ng-l)  exp 


( 


y-£ 

E 


(A10) 


Now,  since  l  > >  E  due  to  the  smallness  of  Ln  relative  to 
-2  A 

r  (i.e..  L  ■  10  cm  whereas  l  =  10cm),  n(0)  ~  1  ;  furthermore 

n(i)  =  n  For  values  n,  >  0.5  ,  it  is  easy  to  see  that  the 

s  s 

above  approximation  is  valid.  Results  shown  in  a  previous 
section  indicate  that  this  is  reasonably  true.  Thus,  diffusion 
is  unimportant  here. 

In  dimensional  form,  Equation  A10  becomes 


n 

n  * 


1 


,n*  \ 

(~ — 1)  exp 


( - 

eFLn 


y-fe 

2/kTQ 


-1 


n*  =  v^/a 


L 


n 


vi 


(All  ) 
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CENTER 


If  we  now  proceed  to  re-examine  Equations  A1  and  A2 ,  assuming 


that  the  magnitude  of  n"  and  q  are  insignificant  but  that  the  E- 
beam  is  operating,  we  have 


eD .  E  , 

1  dn 

kT. 

1 


dy 


he  +  s 


eD  E  , 
e  dn 

dy 


he  +  S 


( A12 ) 


(A13) 


Clearly,  the  above  set  of  equations  has  only  a  trivial,  non¬ 
physical  solution.  We  see  again  that  a  one  dimensional  forulation 
is  not  feasible. 

It  was  mentioned  in  Section  II-l  that  nc  ^  0  was  a  require¬ 
ment  for  stability  in  the  electrode  regions.  We  now  briefly  dis¬ 
cuss  what  this  implies  in  the  presence  of  the  E-beam.  Recall 


ne  =  vi  -  an  +  S 


(A14) 


Let 


n  =  S  /a 
o  ' 


So  that 
Solving  for  n,  a  f  0  , 


2  2 

ne  =  a(n*n  +  nQ  -  n  )  3>  0 


(A15) 


n  < 


in*  +yr~  +  no 2 


( A16 ) 


Now  n*  is  a  function  of  E/M  as  seen  from  Figure  1  and 


typically  n*  is  much  less  than  nQ  away  from  the  electrode 


According  to  Equation  A16,  potential  difficulties  can 


regions . 

be  encountered  in  regions  where  S  is  sufficiently  small  so 

that  as  E/N  and  n*  decrease  away  from  the  electrode  (i.e.  ,  as 

the  discharge  becomes  non-self-sustaining),  n  <  0  .  This  is 

because  n  can  exceed  n  locally  at  the  fringes  of  the  E-beam 

o 

Whether  or  not  such  a  thing  can  trigger  arcing  remains  to  be 
est  abl ished . 
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APPENDIX  B 


MORE  ON  PROBLEM  DIMENSIONALITY 

The  arguments  introduced  in  Section  II. 3  relating  to  the 
dimensionality  Equations  11-13  are  based  on  physical  grounds. 
Here,  we  present  two  more  mathematical  proofs  which,  while  not 
sufficient,  show  clearly  the  inappropriateness  of  the  resulting 
solutions.  It  is  possible  to  combine  Equations  16-18  into  the 
single  equation  for  the  electric  field  seen  below 


kVo,,d3E  kVo  d  r  dE  kTo  +  dE^  co  d  ,„2, 

e'E  dy 3  ’  2e“E2  dy  L  dy  J  eE^  k  2kT  Hv'h  ' 


'dy'  2kTQ  dy 


K  =  0 


where  K+  =  ^/eD.  +  Je/eDe  and  K  = 
From  Reference  12  5  =  n.  ~ne=e^  dE/dy 


L  ,  r 

a  =  n,  +  n  =  +  K  ) 

i  e  eE  dy 


Equation  B1  is  a  third  order,  ordinary  non-linear 
differential  equation.  We  may  nondimensionalize  it  with 


E  =  _ 1 

o  eL 


e  kT 
o  o 


E  «--*•  E/E 


y  ^y/L 


From  a  dimensional  analysis  of  Equation  B1  ,  it  is 
found  that  the  above  substitutions  will  result  in  a  form  of 
the  equation  free  of  coefficients.  Equation  B1  thus  becomes, 


substituting  primes  for  the  derivatives 


E '  '  '  -  ^(E'2)’  -  |  E'  -  |(E2)'  +  EK  =  0  (B3 ) 

Where,  in  the  above,  E,  K,  and  y  are  nondimensional .  We  can 
look  for  solutions  to  this  equation  which  monotonically 
decrease  from  the  electrode  to  a  small  value  at  the  undisturbed 
plasma . 

A  solution  which  is  both  very  simple  and  reasonable  for 
the  anode  is 

E(y)  =  2(y  +  a)-1  (B4) 

a  =  2/Ea 


This  solution  satisfies  Equation  B3  for  the  special  case 


K+ 

=  -2K 

(B5) 

or 

h 

1  Di  . 

-  3  Be  Je 

(B6 ) 

Moreover , 

E' 

=  6  =  ni  -  ne 
=  -2 (y  +  a)-2 

(B7) 

is  also  reasonable  for  the 

anode  since  the 

space 

charge 

is 

negative  and  rapidly 

vanishes  as  y  grows. 

Now 

the  sum 

of 

the  charge  densities  is  found  from 

°  =  ni+ne  =  I[K++n|]  (B8) 

=  [K++4( y+a)-2]  (B9) 


•  -  TufajaWfirfatiit. 
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But  this  result  is  unreasonable  since  a  must  be  zero  (or 
near  zero)  at  the  electrode  and  increase  monotonically  to  an 
asymptotic  value  at  the  undisturbed  plasma.  The  result  shown 
in  Equation  B9  is  clearly  not  a  physically  acceptable  one. 

The  solution  for  the  cathode  is  similarly  obtained  with 


E  =  2(a-yT1,  y<a  (BIO) 

and  K+  =  2K~  (Bll) 

Please  note  that  the  cathode  is  to  the  right  of  the  anode,  as 
indicated  in  Figure  2  of  the  text.  In  the  above  solution, 
both  E  and  6  appear  to  be  reasonable  (the  space  charge 
being  positive  at  the  cathode).  But,  again,  a  does  not 
follow  the  physically  acceptable  pattern. 

A  more  complete  solution  to  Equation  BE  may  be  obtained 
if  we  note  that  since  K~"  is  zero  in  the  undisturbed  plasma, 
it  follows  that  it  must  be  zero  everywhere  as  long  as  we  insist 
on  a  one-dimensional,  constant  property  solution  to  Equations 
16-18.  We  have,  multiplying  p.3  by  E, 

2 

EE'''  -i(E'2)’  -  K+E '  -  |_(E2)'  =  0  #  K+  =  2je/eDc  (B12) 


with  E  ->  E  , 

oo 

E’  -*•  0,  E '  ' 

-*•  0  as  y 

->  the 

undisturbed  plasma 

Now 

EE  '  '  '  = 

(EE'  '  )  ' 

,2), 

E2  2  , 

1,4, 

and 

|-(E  )' 

=  4(E  } 

subst  ituting  into  T>12 

i 

[EE'  '  -  (E'  2  -  K+E  -  i-E1  ]  =  0 
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Int egrat ing 

E"  -  (E')2  -  K+E  -  Je4  =  Cj  (B13 ) 

We  evaluate  C-^  at  y  -*■  the  undisturbed  plasma 


thus 


+  14 

C  =  -K  E  -  4e 
2  00  4  00 


EE"  -  (E‘)2  -  K+(E-Eco)  -  i(E4-E J)  =  0  (B14) 


After  some  algebraic  manipulation  we  find  expressions  for 
and  E"  , 


E' 


[K+(-2E  +  Em  +  §-) 

00 


+ 


2E0O2E2 ) 


(B15 ) 


E"  =  i[2K+(-l  +  |  )  +  (E3  -  Eoo2E ) ]  (B16) 

CO 

E'"  =  \  [  2K+/Eoo  +  (3E2  -  Em2 )  ]  E'  .  as  E+E^ ,  E'-O  and 
E "  '  -+  0  as  well . 

It  is  possible  to  solve  for  a  at  this  point 


Furthermore , 


+ 


K+  =  a  so  that 

00  00 


a 


+ 


) 


(B17) 


(B18  ) 


Clearly,  since  E  _>  EM>  o  _>  oB  which  is  an  improper  be¬ 
havior  for  the  sum  of  the  charge  densities.  Again,  while  6 
and  E  seem  to  behave  reasonably,  the  one-dimensional  solu¬ 
tion  is  incorrect  for  a  . 
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APPENDIX  C 


COMPUTED  PROGRAM 


This  appendix  includes  the  program  listing  with  some  pre¬ 
liminary  comments. 

As  previously  discussed  the  matrix  form  of  the  iteration 
scheme  is:  Rth  Coefficients 


Ma  t  r  i  x 


1  rA,ke 


There  are  three  equations  to  be  solved  at  each  (  i  , , j  ) 

grid  point  corresponding  to  ,  Fg  and  (equations  27,  28  and  29) 
The  solution  vector  is  interlaced  as: 

(  ,  j),  Ane(i,  j),Ani(i,  j),  A<J>(i,j+l)... _ J 


This  method  was  chosen  to  keep  the  values  of  neighboring  points 
in  the  two  dimensional  grid  mesh  near  each  other  in  the  solution 


vector.  Since  a  grid  point  value  Fy 


can  be  expressed  as  a 


linear  combination  of  nearest  neighbor  4>(  i  +1 ,  j+l )  ,  ne(i  +  l,j+l) 

or n^(i+l  ) , j  +  1 )  this  would  produce  a  banded  coefficient  matrix 
A.  of  width  46  and  length  273  (for  a  7x13  grid  mesh)  when  eco¬ 
nomically  stored.  The  International  Mathematics  and  Statistics 
Library  (IMSL)  subroutine.  LEQT1B ,  is  suited  for  solving  equa¬ 
tions  of  this  form. 

The  Sheath  program  calculates  the  coefficients  of  the 
matrix  A  and  vector  F  (refered  to  as  vector  C  in  the  computer 
program).  Thus  C(3n-2)  is  the  F^  value  at  grid  mesh  (i,j), 


C(3n^l)  is  the  value  of  the  Fg  value  at  grid  mesh  i,j  et 
(where  n=i+j-l).  The  solution  to  the  matrix  equation 
[A<J>(i,j),  etc.]  is  assigned  to  vector  C  upon  return  from 
LEQT1B  subroutine. 

The  solution  set  4>^ ,  n& ,  and  n^  are  updated  per 
Equation  31.  The  electron  energy  term  8  .  is  programmed  as 
empirical  function  of  E/N  at  grid  point  (i,j)  and  is  cal¬ 
culated  at  the  completion  of  each  iteration.  The  iteration 
sequence  is  repeated  as  many  times  as  required  to  converge 
to  a  solution. 

The  program  listing  follows. 
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